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Abstract. Tomographic and 3-D analyses for extended, emission-line objects are applied to long-slit ESO NTT -|- 
EMMI high-resolution spectra of the intriguing planetary nebula NGC 7009, covered at twelve position angles. We 
derive the gas expansion law, the diagnostics and ionic radial profiles, the distance and the central star parameters, 
the nebular photo-ionization model and the spatial recovery of the plasma structure and evolution. The Saturn 
Nebula (distance~1.4 kpc, age~6000 yr, ionized mass~0.18 Mq) consists of several interconnected components, 
characterized by different morphology, physical conditions, excitation and kinematics. We identify four "large- 
scale", medium-to-high excitation sub-systems (the internal shell, the main shell, the outer shell and the halo), 
and as many "small-scale" ones: the caps (strings of low-excitation knots within the outer shell), the ansae (polar, 
low-excitation, likely shocked layers), the streams (high- excitation polar regions connecting the main shell with 
the ansae), and an equatorial, medium-to-low excitation pseudo-ring within the outer shell. The internal shell, 
the main shell, the streams and the ansae expand at Vcxp — 4.0xR" km s~^, the outer shell, the caps and the 
equatorial pseudo-ring at V^cxp — 3.15xR" km s~^, and the halo at V^cxp —10 km s~^. We compare the radial 
distribution of the physical conditions and the line fluxes observed in the eight sub-systems with the theoretical 
profiles coming from the photo-ionization code CLOUDY, inferring that all the spectral characteristics of NGC 
7009 are explainable in terms of photo-ionization by the central star, a hot (logT«~4.95) and luminous (log 
L*/L0~3.7O) 0.60-0.61 M0 post-AGB star in the hydrogen-shell nuclear burning phase. The 3-D shaping of the 
Saturn Nebula is discussed within an evolutionary scenario dominated by photo-ionization and supported by the 
fast stellar wind: it begins with the superwind ejection (first isotropic, then polar deficient), passes through the 
neutral, transition phase (lasting ~ 3000 yr), the ionization start (occurred ~2000 yr ago), and the full ionization 
of the main shell (~1000 yr ago), at last reaching the present days: the whole nebula is optically thin to the UV 
stellar fiux, except the caps (mean latitude condensations in the outer shell, shadowed by the main shell) and the 
ansae (supersonic ionization fronts along the major axis). 

Key words, planetary nebulae: individual: NGC 7009- ISM: kinematics and dynamics 



1. Introduction 

NGC 7009 (PNG 037.7-34.5, Acker et al. 1992) is a 
fairly bright, well-studied planetary nebula (PN) exhibit- 

* Based on observations made with: ESO Telescopes at 
the La Silla Observatories (program ID 65.1-0524), and the 
NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Institute. Observing programs: 
GO 6117 (P.I. Bruce Balick), GO 6119 (P.I. Howard Bond) and 
GO 8390 (P.I. Arsen Hajian). STScI is operated by the asso- 
ciation of Universities for Research in Astronomy, Inc. under 
the NASA contract NAS 5-26555. We extensively apply the 
photo-ionization code CLOUDY, developed at the Institute of 
Astronomy of the Cambridge University (Ferland et al. 1998) . 
Send offprint requests to: F. Sabbadin, sabbadin@pd.astro.it 



ing large stratification effects. It is also named the "Saturn 
Nebula" because of the peculiar optical appearance, as 
seen in Fig. 1, where the principal morphological sub- 
systems are indicated. The high-excitation He II emission 
forms a faint, elongated, irregular disc with an oval ring 
better defined out of the apparent major axis. In Ha the 
deformed, hollow main shell (sometimes called "the rim" ; 
Balick et al. 1994, 1998) presents two faint streams spread- 
ing out along the major axis; it is surrounded by an outer 
shell and an even weaker, diffuse halo. The [N II] emission 
(bottom panel in Fig. 1; the gaps in the WFPC2 detectors 
are clearly visible) is concentrated in the external ansae 
along the major axis, and in the caps (i.e. two series of 
condensations not aligned with the major axis). 
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Fig. 1. WFPC2 (HST; programme GO 6117; P.I. Bruce 
Balick) appearance of NGC 7009 (logarithmic scale) in 
the high-excitation line A4686 A of He II (top panel), in 
Ha (which is emitted by the whole ionized nebula; central 
panel), and in the low-excitation line A6584 A of [N II] 
(bottom panel; the image is spread in the four WFPC2 
detectors). North is up and East to the left. The main 
morphological sub-systems are indicated. 

In short: the main and outer shells, the halo and 
the streams connecting the main shell with the ansae 
are medium-to-high excitation regions, whereas the caps 
and the ansae present enhanced low-to- very low-excitation 
emissions. Spectacular HST multi-color reproductions of 
the Saturn nebula are in Balick et al. (1998) and Hajian 
& Terzian at http: / / ad.usno.navy| mil/pne/ gallery. html. 

Based on high-dispersion spectroscopy, Weedman 
(1968) argues that the main shell of NGC 7009 consists 
of a thin prolate spheroid (an ellipse rotating around the 
major axis) seen perpendicular to the major axis. 

Fabry-Perot interferometry and narrow-band imaging 
by Reay & Atherton (1985) indicate that the caps and the 
ansae constitute a double pair of condensations symmetri- 
cally displaced in velocity about the central star. The caps 
expand at ~38 km s^^ (the western condensation is ap- 
proaching and the eastern one receding); the line through 
the caps has an angle of 51° with the line of sight. The 
ansae move at ~60 km s^^ (the eastern one is approach- 
ing and the western one receding); the line of sight forms 
an angle of 84° with the line through the ansae, i.e. they 
almost lie in the plane of the sky. 

Ground-based spectroscopy and HST multi- 
wavelength imagery by Balick et al. (1994, 1998) 
indicate that the caps and the ansae are genuine N- 
enriched FLIERs (fast, low- ionization emission regions; 
according to the definition by Balick et al. 1993), each cap 
consisting of a complex network of small knots. Although 
the "hammerhead" appearance at the outer edge of the 
ansae is typical of a shock front, both pairs of FLIERs 
exhibit a decreasing ionization with distance from the 
star, as expected of a photo-ionized gas. 

Also according to Hyung & AUer (1995a, b) the ob- 
served spectrum of NGC 7009 is fairly well fitted by a 
photo-ionization model. They analyze two regions in the 
bright ring on the major and minor axes (dominated by 
the low-excitation and the high-excitation lines, respec- 
tively), obtain that the caps are N-rich knots and the 
ansae, very likely, do not lie along an axis almost per- 
pendicular to the observer, and suggest a nebular model 
consisting of a dense toroid ring (or quasi-fiat shell) + 
broad conical shells. 

In a very recent paper Gongalves et al. (2003) present 
low-resolution spectra along the major axis of NGC 7009, 
and discuss the physical excitation and chemical proper- 
ties of the main shell, the caps, the streams and the ansae, 
obtaining that: 

- all regions are mainly radiatively excited, 



- there are no clear abundance changes across the neb- 
ula, but only marginal evidence for a modest overabun- 
dance of nitrogen in the ansae, 

- none of the available theoretical models is able to 
account for the observed characteristics of the ansae. 

Gongalves et al. (2003) also stress the need for an ac- 
curate kinematical study of the nebula. 

Deep [O III] imaging by Moreno-Corral et al. (1998) re- 
veals the presence of a large (r~130"), faint, complex halo 
around NGC 7009, characterized by a series of fractured 
filaments, semi-circular envelopes and irregular condensa- 
tions. 

Our nebula is an extended X-ray emitting source 
(Guerrero et al. 2002): its central cavity consists of low- 
density (a few tens H-atom cm~^), hot (T~1.8xl0^ K), 
shocked fast stellar wind. In fact, the exciting star is losing 
mass at high velocity, Vcdgc— 2700 km s~-^ (Cerruti-Sola & 
Perinotto 1985) and at an uncertain rate: 3.0x10"^" M0 
yr-i (Cerruti-Sola & Perinotto 1985), 1.0x10"® Mq yr'^ 
(Bombeck et al. 1986), 2.8x10"^ Mq yr"! (Cerruti-Sola 
& Perinotto 1989), 2.1x10"^ Mq yr^^ (Hutsemekers & 
Surdej 1989), and 2.0x10"^ Mg yr"! (Tinkler & Lamers 
2002). 

In spite of the great deal of attention, both theoret- 
ical and observational, and of the rich bibliography on 
NGC 7009, many (even fundamental) topics remain un- 
certain, e.g.: (a) the nebular distance (between 0.6 and 
2.5 kpc), and the central star mass (0.54 Mq <M* <0.70 
M0) and luminosity (3.0<log L^/L© <4.2); (b) the ex- 
citation processes (shocks or photo-ionization) and the 
physical mechanisms (ionization-front instabilities, stel- 
lar or nebular bullets, photo-evaporation of neutral knots 
by a diluted, fast stellar wind, Oort-Spitzer rockets, i.e. 
photo-evaporation of neutral knots in the absence of a fast 
stellar wind; for details, see Balick et al. 1998) forming 
and shaping the caps and the ansae; (c) the discrepancy 
in the chemical abundances obtained from recombination 
and collisionally excited lines (Hyung & AUer 1995b, Liu 
et al. 1995, Luo et al. 2001); (d) the overall, de-projected 
spatial structure of the ionized (and neutral) gas making 
up the Saturn Nebula. 

To deepen our insights in the kinematics, physical con- 
ditions, ionic and spatial structure, and evolutionary sta- 
tus of this fascinating object, NGC 7009 was included in 
the sample of proto-PNe and PNe observed with ESO 
NTT-f EMMI (spectral range AA3900-7900 A, spectral res- 
olution A/AA=R=60 000, spatial resolution S~1.0") and 
the Telescopic Nazionale Galileo (TNG)+SARG (spectral 
range AA4600-8000 A, R=115 000, S~0.7"). The echel- 
lograms were reduced and analyzed according to the 3-D 
methodology so far applied to NGC 6565 (Turatto et al. 
2002) and NGC 6818 (Benetti et al. 2003). 

The results are contained in this paper, whose struc- 
ture is as follows: Sect. 2 describes the observational 
procedure, the spectroscopic material and the reduction 
method. Sect. 3 disentangles the kinematics of the dif- 
ferent sub-systems. Sect. 4 studies the radial profile of 
the physical conditions (electron temperature and elec- 
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tron density) in different directions, in Sect. 5 we per- 
form a critical analysis of the nebular distance, size, mass 
and age, in Sect. 6 the central star parameters arc given, 
Sect. 7 combines the observed ionic profiles and the photo- 
ionization model (CLOUDY), Sect. 8 describes the 3-D 
spatio kincmatical structure of NGC 7009, and Sect. 9 
contains the general discussion, sketches the shaping of 
the whole nebula and the different sub-systems, and draws 
the conclusions. 

2. Observations and reductions 

To obtain an adequate spatial and spectral coverage of 
this extended and complex nebula, NGC 7009 was ob- 
served with ESQ NTT + EMMI (cchcUc mode; grating 
#14, grism #3) at twelve equally spaced position angles 
(PA), including PA=79°, the apparent major axis con- 
taining the streams and the ansae. In all cases the spec- 
trograph slit (1.0" wide and 60" long) was centered on the 
exciting star, whose continuum is used as position marker. 
During the observations the sky was photometric and the 
seeing fluctuated between 0.50" and 0.75". 

The echellograms (exposure times 360s to 600s) 
cover the spectral range AA3900-7900 A with resolution 
R~60 000, and provide the spatio-kinematical structure of 
the main ionic species ([O I], ionization potential (IP)=0 
eV at A6300 A to [Ar V], IP=59.8 cV at A7005 A) within 
the nebular slices covered by the spectrograph slit. 

Bias, zero-order flat field and distortion corrections, 
and wavelength and flux calibrations were performed ac- 
cording to the straightforward procedure described by 
Turatto et al. (2002). For educational purposes in Figs. 
2 and 3 we present the spectral structure of NGC 7009 
in each of the twelve observed PA at low, mean and high- 
excitation (A6584 A of [N II], A5007 A of [O III] and A4686 
A of He II, respectively), and at Ha (A6563 A of H I), 
which is emitted through the whole of the ionized nebula, 
independent of the plasma excitation conditions. 

These pictures highlight the complex ionization struc- 
ture of the Saturn Nebula: 

- the N"*" emission is faint everywhere, except in the caps 

and the ansae. The same is observed at even lower ex- 
citation, like S+ at A6717 A and A6731 A, and 0° at 
A6300 A. Note that, due to the large optical extent of 
the nebula along the major axis (PA=79°), in Fig. 2 
(bottom panel) the weak [S HI] line at A6312 A (echelle 
order 97) partially overlaps the A6584 A [N II] emission 
(echelle order 93) of the W-SW ansa. Hereafter the net 
contribution of [N II] at this PA is obtained thanks to 
the similarity of the ionization potential of S~^~^ and 
Ar++ (23.4 eV and 27.6 eV, respectively): we subtract 
the whole spectral image; of [Ar HI] at A7135 A (prop- 
erly scaled in flux) from that of [S HI] at A6312 A, as in 
NGC 7009 F(A7135 A)/F(A6312 A)~10.0 (Kingsburgh 
& Barlow 1994, Hyung & AUer 1995a, b, this paper); 

- Ha is blurred by thermal motions, fine structure and 

expansion velocity gradient. Although this masks the 
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detailed structure of the ionized gas, the H+ distribu- 
tion mimics that of 0++; 

- A5007 A of O"*""*", the strongest line in the optical region, 

represents the overall characteristics of the ionized gas 
(except the innermost regions, where 0+^ and higher- 
ionization species prevail, and the outer edge of the 
caps and the ansae, rich in 0+ and 0°); 

- the He"*" A4686 A line marks the highest-excitation neb- 

ular regions. The weak, blue-shifted, diffuse emission 
present at all PA is due to the tail in the thirteen fine- 
structure components of this recombination line (also 
broadened by thermal motions). 

From Figs. 2 and 3 we infer that the main shell of 
NGC 7009 is a thin tri-axial ellipsoid, broken (or frac- 
tured) along the major axis (projected in PA=79°), where 
the streams and the ansae appear. Note the line tilt along 
the apparent minor axis (PA=169°), suggesting that we 
are not aligned with either the intermediate or the minor 
axes of the ellipsoid, and that the line of the nodes lies at 
PA~130° 310°. The outer shell is an attached, complex, 
broad envelope completely filled with matter, as suggested 
by the absence of splitting at the rest wavelength of the 
[O HI] line (out of the main shell) . Kinematically the caps 
belong to the outer shell. 

Note that the main shell is red-shifted (by a few 
km s~^) with respect to the outer shell (-1- caps) and 
the halo; when considering the baricentric radial veloc- 
ity of the ansae and the caps, we have: T^ad (outer shell) 

- Frad (ansae) ~ Vrad(caps) V;ad(halo) ~ Vrad(main 
shell) - 3(±1) km s'^ The streams at PA=79° present 
hybrid kinematics: the internal region is co- moving with 
the main shell, the external one with the ansae (and the 
halo). Also note the variable FWHM along the streams, 
and the hooked appearance of the ansae (with both "tips" 
blue-shifted!). 

Our preliminary, qualitative analysis stresses the com- 
plexity of NGC 7009. This appears to be a common char- 
acteristic of all PNe covered at adequate spatial and spec- 
tral resolutions, and questions the validity of the spatio- 
kinematical studies based on unsuitable observational ma- 
terial or, even worse, on the assumption of spherical sym- 
metry for the emitting gas. Moreover, Figs. 2 and 3 
can only provide a rough picture of the large amount of 
kinematical and physical information stored in the echel- 
lograms. Let's start our thorough investigation with the 
gas kinematics. 

3. The gas spatio-kinematics 

From Fig. 4, the spectrum of an elementary volume within 
the slice of nebula selected by the spectrograph slit is char- 
acterized by: 

Vrad = Vneh " Kxp X COS a (1) 

and 

r = \R X sina|, (2) 
where: 
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Fig. 2. Detailed structure (on a logarithmic scale) of representative emission at the twelve observed PA of NGC 7009. 
The original fluxes are multiplied by the factor given in parentheses, to render each emission comparable with A5007 
A of [O III]. The blue-shifted gas is to the left. The slit orientation at each PA is indicated in the He II frame. At 
PA=79° (apparent major axis) the weak [S III] line at A6312 A (echelle order 97) partially overlaps the A6584 A [N 
II] emission (echelle order 93) of the W-SW ansa; see the text for details. 



Fig. 3. Same as Fig. 2. 

Vnch is the systemic radial velocity of the whole nebula; 
14xp is the expansion velocity; for PNe it is generally 

written in the form T4xp=cixR+C2 (Wilson 1950, 

Weedman 1968); 
R and r are the true and the apparent distances from the 

central star, respectively; 
a is the angle between the line of sight and the central 

star-elementary volume direction. 

Following Turatto et al. (2002, and references therein), for 
each spectral image of NGC 7009 we select: 

- the "central star pixel line" (cspl), representing the neb- 

ular gas projected at the apparent position of the star, 
whose motion is purely radial. The cspl is the same 
at all PA; in Eqs. (1) and (2) it corresponds to a=0° 
and 180°; i.e. 2yexp=Vrad(a=180°) - Frad(a=0°), and 
rcspi=0; 

- the "zero velocity pixel column" (zvpc), which is unaf- 

fected by the expansion velocity field, and gives the 
spatial distribution of the ionized gas moving per- 
pendicular to the line of sight (valid for a=90° and 
270°; i.e. K-ad=Keb, and 2r^vpc=2Rzvpc =R(a=90°) + 
R(a=270°)). 

We will derive the overall spatio-kinematical properties of 
the Saturn Nebula by assembling the kinematical results 
in a single radial direction (the cspl), and the spatial re- 
sults in twelve tangential directions (the zvpc at each PA). 
Given the nebular complexity, the different sub-systems 
will be analyzed separately. 

3.1. Kinematics of the main shell (+ streams + ansae) 

The peak separations in the cspl, 2Vcxp, are contained in 
Table 1 (last column), where the ions arc put in order 
of increasing IP ([O I] at A6300 A is absent, since the 
main shell is a medium-to-high excitation region). Typical 
errors are 1.0 km s~^ for the strongest forbidden emissions 
(like A5007 A of [O III]) to 2.0 km s'^ for the faintest ones 
(in particular, A5517 A and A5537 A of [CI III], A6578 
A of C II and A7005 A of [Ar V]). The errors for the 
recombination lines are: 2.0 km s-i for A6563 A of H I, 
and 1.5 km s'^ for A5876 A of He I and A4686 A of He II. 
Table 1 also contains the kinematical results taken from 
the literature. 

Concerning the zvpc in the main shell, the intensity 
peak separations, 2rzvpc, in the different emissions at the 
twelve PA of NGC 7009 are Usted in Table 2. Please note 
the different "relative" spectral and spatial resolutions 



Fig. 5. The complex kinematics of NGC 7009 at PA= 
169° (apparent minor axis). Empty circles= zvpc vs cspl 
in the main shell, providing V|3xp(main shell) (km s~^)= 
4.0(±0.3) xR"; dots= extent of both the zvpc and the cspl 
at 10% maximum intensity; this is valid for the outer shell 
and gives Kxp (outer shell) (km s-i)= 3.15(±0.3) xR"; 
moreover. Vcxp(halo)~10 km s^^. Long-dashed line= ex- 
pansion law from Weedman (1968); continuous line= ex- 
pansion law adopted in this paper. 

of our echellograms: the former is given by yexp/AV~4 
(AV=spectral resolution), and the latter by r/Ar ~8-16 
(r=apparent radius, Ar=seeing); this means that the spa- 
tial information of NGC 7009 is more accurate than the 
kinematical information. 

According to Benetti et al. (2003, and references 
therein), the next step is the identification of the PA 
at which R^vpc— Rcspi- Let us assume for the main shell 
of NGC 7009 the most general spatial structure, i. e. a 
tri-axial ellipsoid. The major axis, projected in PA=79°, 
is almost perpendicular to the line of sight, since the 
emissions appear nearly un-tilted. The intermediate and 
the minor axes lie close to PA~169°; actually, the line- 
tilt observed here suggests that we are not aligned with 
either axis (see Fig. 4). We conclude that R(minor 
axis)<Rzvpc— Rcspi<R(intermediate axis) along the ap- 
parent minor axis of the nebula, i.e. at PA=169°. 

The cspl-zvpc trend for the main shell of NGC 7009 
at PA=169°, shown in Fig. 5, follows the Wilson law: 
the high-excitation particles expand more slowly than the 
low-excitation ones. Moreover, there is a positive correla- 
tion between the expansion velocity and the size of the 
monochromatic image (Wilson 1950). Note that, in spite 
of the wide range of IP observed (10.4 eV of S"*" to 59.8 
eV of Ar+*), the range of both 2rzvpc and 2V^xp is quite 
limited (8.3 to 10.5 arcsec and 33.0 to 41.9 km s^^, respec- 
tively). When added to the distribution of the data in Fig. 
5, this suggests that NGC 7009 at PA=169° is optically 
thick in Ar+^ and Hc++, nearly thick in N^\ Ne^'^ and 
Ar+'^, and thin in the other ionic species, and/or that the 
radial density profile of the main nebula is quite sharp. To 
gain insight into the question, a kinematical study at even 
higher ionization states appears indispensable, e.g. in the 
forbidden line of Ne+^ at A3425 A (which, unfortunately, 
is outside our spectral range). 

The r^vpc vs IP trend being the same at all observed 
PA (from Table 2), we can assess that the expansion ve- 
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Fig. 4. Nebular slice vs. spectral image connection for a typical PN. In this example we use the spectral image of 
the [O III] A5007 A line of NGC 7009 (brightest layers of the main shell at PA=169°; see Fig. 3, bottom panel). The 
"central star pixel line" (cspl) and the "zero velocity pixel column" (zvpc) of the spectral image, as discussed in the 
text, are indicated. 

Table 1. Peak separation in the cspl of NGC 7009 (main shell) 



Ion IP (eV) A (A) 2Kxp (km/s) 

(1) (2) (3) (4) (5) (6) (7) this paper 



[S II] 


10.4 


6717-6731 




- 




- 


- 


41.9 


[0 II] 


13.6 


3727-7319 


40.9 




40.0 






41.5: 


H I 


13.6 


6563 


42.1 


44.0 


35.6 32.4 


36.8 


34.5 


40.0 


[Nil] 


14.5 


6584 












41.6 


[S III] 


23.4 


6312 












41.4 


[CI III] 


23.8 


5517-5537 












41.0: 


C II 


24.4 


6578 












40.0: 


He I 


24.6 


5876 


44.1 










40.5 


[Ar III] 


27.6 


7135 












40.8 


[0 III] 


35.1 


5007 


41.1 




41.2 






40.8 


[Ar IV] 


40.7 


4711-4740 


39.7 










39.7 


[Ne III] 


41.0 


3869-3967 


38.7 










40.0 


N III 


47.4 


4640 


41.1 










40.0: 


He II 


54.4 


4686 


35.8 




34.6 






36.6 


HI 


54.9 


3444-3770 


32.9 








30.5 




[Ar V] 


59.8 


7005 












33.0: 


(1)= 


Wilson (1950) 




(5)= 


Cristiani et al. (1989) 








(2) = 


Mcndcz ct al. (1988) 




(6)= 


Bianchi (1992) 









(3) = Meatheringham et al. (1988) (7)= Liu & Danziger (1993) 

(4) = Bianchi et al. (1989) 

Table 2. Peak separation in the zvpc at the twelve PA of NGC 7009 (main shell) 



Ion 2rzvpc (arcsec) 





PA=04° 


19° 


34° 


49° 


64° 


79° 


94° 


109° 


124° 


139° 


154° 


169° 


[S II] 


12.1 


13.4 


14.4 


18.5 


23.5 


24.8 


19.4 


13.6 


11.4 


10.6 


10.6 


10.5 


[0 II] 


12.0 


13.2 


14.2 


18.0: 








13.4 


11.2 


10.4 


10.4 


10.3 


HI 


11.2 


12.8 


13.6 


17.0 


22.5 


24.5 


18.5 


12.7 


10.7 


10.0 


10.0 


10.2 


[N II] 


12.1 


13.3 


14.3 


18.3 


23.5 


24.7 


19.3 


13.4 


11.3 


10.5 


10.4 


10.4 


[S HI] 


11.9 


13.1 


14.3 


18.1 


23.3 


24.5 


19.0 


13.3 


11.2 


10.3 


10.3 


10.4 


[CI HI] 


12.0: 


13.0: 


13.8: 


18.2: 








13.0: 






10.3: 


10.4: 


C II 


11.5: 


13.1 


14.1: 


18.0: 








12.5: 


10.8: 


10.0: 




10.4: 


He I 


12.0 


13.0 


14.1 


18.0 


23.5 


24.5 


18.7 


13.2 


11.2 


10.3 


10.3 


10.4 


[Ar HI] 


12.0 


13.2 


14.0 


18.1 


23.4 


24.5 


18.5 


13.2 


11.2 


10.3 


10.2 


10.3 


[0 HI] 


12.0 


13.2 


14.1 


18.0 


23.3 


24.5 


18.3 


13.0 


11.2 


10.3 


10.3 


10.2 


[Ar IV] 


11.3 
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locity in the main shell of NGC 7009 is proportional to 
the distance from the central star through the relation: 

VexpCmain shell)(km s"^) = 4.0(±0.3) x R" . (3) 

Weedman (1968) has obtained a different law: T4xp(km 
s~^)= 7.1 X {R" — 6.9) (also shown in Fig. 5), which ap- 
pears questionable, since it is based on: (I) spectra at a 
single PA (along the apparent major axis); (II) the "a pri- 
ori" assumption that the nebula is a prolate spheroid with 
a=27.7" and a/b=1.5, seen perpendicular to the major 
axis. 

Although weakened by projection effects, Eq. (3) is 

extended to the streams and the ansae, for reasons of 
continuity (Figs. 2 and 3). Both these sub-systems ap- 
pear aligned to the major axis of the central ellipsoid, and 
highly inclined with respect to the line of sight (the result- 
ing supersonic velocity of the ansae, l^xp (ansae) ~ 100 km 
s^^, is quite uncertain). Note the large, variable FWHM 
along the streams, and the hammerhead appearance of the 
ansae, suggestive of local turbulent motions (i.e. interac- 
tion with the slower outer shell and halo). 

3.2. Kinematics of tine outer shell (+ caps) 

The foregoing analysis cannot be applied to the outer 
shell, the cspl of NGC 7009 being dominated by the bright 
main shell. In order to derive the spatio-kinematical prop- 
erties of the outer shell, we measure the extent of both 
the cspl and the zvpc at 10% maximum intensity (this 
cut avoids the contribution of the halo to the zvpc) in the 
strongest lines (A5007 A of [O III], A5876 A of Ho I, A6584 
A of [N II] and A7135 A of [Ar III]) at PA=169° (the ap- 
parent minor axis). The results arc graphically shown in 
Fig. 5. When extended to the other PA, and combined 
with the whole emission structure of the outer shell, they 
furnish the following expansion law: 

14xp(outer shell)(km s"^) = 3.15(±0.3) x R" , (4) 

valid for the caps too. 

3.3. Kinematics of the halo 

Only at PA=64° and 94° does the faint [O III] A5007 A 
emission of the halo present a two-peaks velocity struc- 
ture, providing yoxp(halo) = 12(±2) km s^^, whereas it 
is not spUt at the other PA, with FWHM~2T/cxp(halo) 
= 20(±5) km s-^ The line split observed at PA=64° and 
94° (i.e. close to the apparent major axis) is connected to 
the extended emissions surrounding the streams and the 
ansae (see Fig. 1), which are reminiscent of bow shocks 
(Balick et al. 1998). 

3.4. Kinematics overview 

The intriguing kinematical properties of the Saturn 
Nebula are presented in Fig. 6, showing the position- 
velocity (P-V) maps, i.e. the complete radial velocity field 



at the twelve observed PA. We have selected He II, [O 
III] and [N II] as markers of the high, medium and low- 
excitation regions, respectively. These P-V maps are rela- 
tive to the systemic heliocentric velocity of the main shell, 
VJ.ad0(main shell)= -46.0(±1.0) km s~^, corresponding to 
^LSR(uiain shell)= -36.3(±1.0) km s~^, and are scaled ac- 
cording to Eq. (3), i.e. they reproduce the tomographic 
maps of the main shell (-|- the streams -|- the ansae) in the 
nebular slices covered by the slit. 

Concerning the outer shell (-|- the caps), the tomo- 
graphic maps in Fig. 6 are: 

(a) compressed along the horizontal axis, since the kine- 
matics of these sub-systems is provided by Eq. (4) ; this 
is particularly evident along and close to the apparent 
minor axis; 

(b) slightly blue-shifted, as V^-ad© (outer shell) ~ 

Viad0(caps) ~ -49 km s"^. 

Note the peculiar behaviour of the ansae: they follow the 
expansion law of the main shell and the streams, whereas 
their baricentric radial velocity is the same as that of the 
outer shell, the caps and the halo. 

Fig. 6 highlights the main limitation of any spatio- 
kinematical study of the Saturn Nebula: the major axis 
(containing the most interesting features, i.e. the streams 
and the ansae) almost lies and moves in the plane of the 
sky, the radial component of the expansion velocity is 
nearly null, and projection effects play the leading role. 
In other words: we are forced to extrapolate the (well- 
defined) kinematical and spatial properties of the matter 
perpendicular to the major axis to the poor kinematics 
and precise spatial profile of the gas along the major axis. 

The P V maps (Fig. 6) stress the variety of kinemat- 
ics in the three "large-scale" sub-systems of the Saturn 
Nebula (the main shell, the outer shell and the halo) and 
in the three "small-scale"" ones (the streams, the caps 
and the ansae) . Different evolutionary scenarios can qual- 
itatively explain the "large-scale" features: (a) a multiple 
(triple) ejection, the halo being the oldest and the main 
shell the youngest; (b) different accelerating mechanisms 
acting on a single, prolongated ejection; (c) a combina- 
tion of (a) and (b). To disentangle these possibilities (and 
their correlation with the "small-scale" structures) is pre- 
mature, requiring more input, and is post-poned to Sects. 
5 to 9. We wish to remark here that the common expansion 
law of the outer shell and the caps suggests that they are 
coeval; very likely, the latter simply represent local con- 
densations within the former (a similar result is obtained 
by Corradi et al. 2000 for the enhanced low-excitation lay- 
ers of IC 2553 and NGC 5882). At present we cannot as- 
sess that the same conclusion is valid for the three other 
co-moving sub-systems of NGC 7009 (i.e. main shell, the 
stream and the ansae). 

Besides the complexity of the kinematics, Fig. 6 high- 
lights the large stratification of the radiation, the Saturn 
Nebula being optically thin in most directions. The ion- 
ization level decreases (i.e. the [N II] emission enhances) 
in: 
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Fig. 7. Overall structure of the [O III] zvpc at the twelve 
observed PA of NGC 7009. The slit orientation is indi- 
cated; moreover: short-dashed line=internal shell; contin- 
uous line=inain shell; long-dashed line=outer shell. 

- two series of knots in the outer shell at PA=64 to 94° 

(i.e. the caps); this appears connected to the shadowing 
by an inner, dense layer (the main shell); 

- the ansae in PA=79°; probably associated to the pres- 

ence of internal, extended, low-density regions along 
the major axis (the streams). 

A final remark concerns the overall structure of the zvpc 
at all PA, providing the matter distribution in the plane 
of the sky passing through the central star. We discard 
Ha as reference emission, because of the large broaden- 
ing, and select A5007 A of [O III] , i.e. the strongest line of 
an ionic species present throughout the whole nebula (in 
first approximation we can put: I(A5007 A)zvpc oc N^"^, 
Ne being the electron density). As shown in Fig. 7, all the 
sub-systems of NGC 7009 are present (except the caps, 
whose contribution in the zvpc is null, due to projection 
effects; see Figs. 2 and 3). Moreover, a further, internal, 
faint and incomplete shell appears close to the major axis 
(PA=49° to PA=109°). As expected, this feature (here- 
after called "the internal shell" ) is also present in the zvpc 
of [Ar III] and, even more strongly. He II, but it is absent 
(or almost absent) in those of [N II] and [S II], i.e. it is 
a high-excitation region. Concerning the expansion law of 
the internal shell, we can only say that it mimics that of 
the main shell. 

In summary, our detailed spatio-kinematical analy- 
sis confirms and extends the morphological complexity 
sketched in Sect. 1, NGC 7009 consisting of a puzzle of 
interconnected sub-systems, as synthesized in Table 3. 

4. The physical conditions 

First of all the observed line intensities must be corrected 
for interstellar absorption according to: 



logg^ = /,c(H/3), 



1(A) 



obs 



(5) 



where f^ is the interstellar extinction coefficient given by 
Seaton (1979). 

The Ha/H/? analysis of the whole spectral image, 
as introduced by Turatto et al. (2002), provides poor 
spatial resolution maps, because of the blurred appear- 
ance of both Ha and H/3. The resulting maps appear 
quite uniform at Ha/H/3=3.20 (±0.10), corresponding to 
c(H/3)=0.15 (±0.05) (for the case B of Baker & Menzel 
1938, Te=10 000 K and log Ne= 3.60; Brocklehurst 1971, 
Aller 1984, Hummer & Storey 1987), with a smooth de- 
cline in the innermost regions (very likely caused by a local 
increase of Te). 

This agrees with most of the c(H/3) estimates reported 
in the literature (0.09-0.27, Perinotto & Benvenuti 1981; 



0.39-0.56, Tylenda et al. 1992; 0.09, Kingsburgh & Barlow 
1994; 0.05-0.28, Bohigas et al. 1994; 0.20, Liu et al. 1995; 
0.09, Hyung & Aller 1995a,b; 0.24, Lame & Pogge 1996; 
0.14, Ciardullo et al. 1999; 0.07, Luo et al. 2001; 0.10, 
Rubin et al. 2002; 0.16, Gongalves et al 2003), indicating 
that NGC 7009 is a PN with little absorption. 

The radial profile of the electron temperature, Te, is 
obtained from diagnostic line ratios (A5007 A/A4363 A of 
[O III] and A6584 A/A5755 A of [N II]), and the radial 
distribution from both line ratios (A6717 A/A6731 A of [S 
II] and A4711 A/A4740 A of [Ar IV]) and the absolute Ha 
flux. We consider the zvpc, corresponding to the gas in 
the plane of the sky, since it is generally independent on 
the expansion velocity field (for details, see Turatto et al. 
2002). In the specific case of NGC 7009 both methods (i.e. 
diagnostic line ratios and Ha flux) suffer some limitations 
due to: 

- the weakness of the [N II] auroral line and the [S II] 

and [Ar IV] doublets. T^N 11], iVo[S II] and A^o[Ar IV] 
can be derived only at the intensity peaks. To get out 
of this trouble we extend the analysis of the diagnostic 
line ratios to the prominent knots of the whole spectral 
images; 

- the large Ha broadening. The deconvolution for instru- 

mental resolution plus thermal motions plus fine struc- 
ture is quite complex, and the resulting F(Ha)z;vpc 
and Ne{H.a) profiles uncertain. Moreover, the diff'er- 
cnt kincmatical properties of the sub-systems must be 
considered, being : 



N^{Ha) oc 



where: 



-0.47 



x( 



ei X rcspi X D 



1/2 



(6) 



- D is the nebular distance; 

- Tcspi is the radius of the cspl in the main shell (this is 

valid for the internal shell, the main shell, the streams 
and the ansae) and in the outer shell (for the outer 
shell and the caps); 

- ei is the "local filling factor" , representing the fraction of 

the local volume actually filled by matter with density 

Thus, according to Benetti et al. (2003), F(Ha)zvpc is also 
obtained from the observed radial ionization structure rel- 
ative to 0++ and the fair assumption that 0/H=constant 
across the nebula (Sect. 7); in the N^. range here consid- 
ered: 



F{Ha) 



zvpc 



F(A5007A) 



zvpc 



CX ^ X /(Te) X icf{0 



(7) 



with ic/(0"'""'")=^qr=ionization correcting factor. 

A comparative analysis gives satisfactory results: the 
iVo(Ha) profiles derived in the two ways differ by less 
than 7%, with the exception of the faint, innermost re- 
gions, where discrepancies as large as 30% are observed. 
At present we are unable to identify the principal factor 
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Fig. 6. Position-velocity maps at the twelve observed PA of NGC 7009 for the high (He II, blue), medium ([O III], 
green) and low ([N II], red) excitation regions, scaled according to the relation \4xp(km s~^)=4.0xR" (vahd for the 
main shell + the streams + the ansae; see the text for details). The orientation of these tomographic maps is as in 
Figs. 2 and 3. 

Table 3. NGC 7009: summary of the kinematical and physical conditions in the sub-systems. 



Sub-system 


Vcxp 


Vrade 




A'e peak 


excitation 


first 




(km s"^) 


(km s"') 


(K) 


(cm-3) 


degree 


citation 


internal shell 


4.0xR"(:) 


-46 (:) 


> 12000 


2000-2500 (:) 


very high 


Sect. 3.4 


main shell 


4.0xR" 


-46 


12000 to 10000 


4000 to 8000 


high 


Sect. 1 


streams 


4.0xR" 


-46 (inner) to -49 (outer) 


10000 (:) 


1000-1500 


high 


Sect. 1 


ansae 


4.0 xR" 


-49 


10000 (:) 


2000-2500 


low to very low 


Sect. 1 


outer shell 


3.15xR" 


-49 


9500-10000 


2000-2500 


high 


Sect. 1 


caps 


3.15xR" 


-49 


9500-10000 


3000-3500 


low to very low 


Sect. 1 


equatorial pseudo-ring 


3.15xR" 


-49 


10000 


2800 


mean to low 


Sect. 8 


halo 


10 


-49 


10000 (:) 


250-300 


high 


Sect. 1 



responsible (inadequate deconvolution? To uncertainties? 
chemical composition gradient?), although the last possi- 
bility is strongly suspected, due to the different nebular 
morphology in IIe+"'" and H+ (Fig. 1); this is discussed 
further in Sect. 7. 

The resulting radial profile rapidly changes with 
PA, as expected of a chaotic object like NGC 7009. Some 

representative examples arc shown in Fig. 8. In order to 
cover all the sub-systems, we take into account: 

- the zvpc at PA=169° for the main shell, the outer shell 

and the halo along the apparent minor axis, 

- the zvpc at PA=79° for the internal shell, the main shell, 

the streams and the ansae along the apparent major 
axis. In the following this direction will be considered 
as the true major axis of the nebula, due to the modest 
line tilt observed at PA=79° (Figs. 2 and 6). 

Concerning the caps, whose contribution to the zvpc is 
null because of perspective effects (Sects. 2 and 3), we 
first de-project the spectral images at PA=94° (through 
Eq. (3) for the main shell, and Eq. (4) for the outer shell 
and the caps), and then consider the direction containing 
the central star and the caps (it forms an angle of 55° with 
the line of sight). 

Let us analyze the physical conditions in the different 
sub-systems (a synthesis is given in Table 3). 



4.1. The internal shell 

We could derive only poor information for this elu- 
sive, high-excitation sub-system: it is a moderate-density, 
A^e(peak)~2000-2500 cm'^, thin region at large Te 
(>12000 K). 



4.2. The main shell 

Tc [O III] presents a well-defined radial profile common at 
all PA: it is > 12 ODD K in the tenuous, innermost regions, 
it rapidly decreases outward down to 10 000 K in the dens- 
est layers, and furtherout it remains more or less constant, 
in quantitative agreement with the results by Rubin et al. 

(2002) , based on HST/WFPC2 imagery. Ground-based 
Te[0 III] determinations (mean value in the brightest re- 
gions) by Kingsburgh & Barlow (1994), Bohigas et al. 
(1994), Hyung & Aller (1995a, b), Liu et al. (1995), Mathis 
et al. (1998), Kwittcr & Henry (1998), Luo et al. (2001) 
and Gongalves et al. (2003) cluster around 10000 K. 

The A'e radial distribution in the main shell shows a 
sharp, bell-shape profile with Nq[S II] peaks up to: 

4000 (±500) cm~^ along the major axis, 
7000 (±500) cm"'' along and close to the minor axis, 
8000 (±500) cm-3 at PA=94°, in the de-projected direc- 
tion containing the caps. 

Previously, Hyung & Aller (1995a, b) obtained iVe[S II] 
peaks of 5000 cm~^ and 5600 cm"-^ along the apparent 
major and minor axes, respectively, and Gongalves et al. 

(2003) 7Ve[S II] =5500-5900 cm'^ along the apparent ma- 
jor axis (cf. Sect. 4.6). 

Ne[Ar IV] is systematically lower than A'^e[S II]. We 
adopt the [Ar IV] electron impact excitation rates given 
by Keenan et al. (1997); when using earlier coUisional rates 
(for example, Aller 1984), 7Ve[Ar IV]~ Ne[S II]. 

In general we have: 7Ve[S H] xef-^ ~ N^iUa) (Aller 
1984, Pottasch 1984, Osterbrock 1989). In Fig. 8 the 
match between A^o (forbidden lines) and Afe(Ha) occurs 
for (rcspixDxei)main shell— 5.0 arcsec kpc, where rcspi=size 
of the main shell in the radial direction, and D=ncbular 
distance (Eq. (6); for details, see Turatto et al. 2002), i.e. 
(Dxei)niain shell =^1-0 kpc (from Table 2). 
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4.3. The streams 

No direct Te determination can be obtained for these 
medium-to-high excitation, low-density (A^e— 1000-1500 
cm"^) regions connecting the main shell with the ansae, 
due to the weakness of the [O III] auroral line. Following 
Rubin et al. (2002) and Gongalves et al. (2003), we adopt 
Te ~ 10000 K. 

4.4. The ansae 

We derive only a rough estimate of re[N II]~ 10 000 K. 
Balick et al. (1994) report T^O III]=11500 K and Te[N 
II]=8100 K, Kwittcr & Henry (1998) To[0 III]=9300 K 
and TelN II]=9100 K, Rubin et al. (2002) T^O III]~9200 
K, and Gongalves et al. (2003) re[0 III] =9300-10 100 K. 
Likely, Tc decreases at the outer edge, where the ion- 
ization drops. Ne[S II] reaches 2000-2500 cm~^ in the 
ansae, i.e. they hardly double the density of the streams, 
in agreement with Balick et al. (1994, 1998), and in par- 
tial agreement with Gongalves et al. (2003), who find Ne[S 
II](streams) ~ N^[S II](ansae)=1300-2000 cm'^. 

4.5. The outer shell 

Our echellograms provide Te[0 III]c:i9500-10000 K and 
Ne[S II] up to 2000-2500 cm^^ Previously, Kwitter 
& Henry (1998) obtained rc[0 III]=9400 K and TefN 
II]=12100 K, and Rubin et al. (2002) T^IO III]=i9200 K. 
Note the presence of a localized knot with Ne[S II] ~ 2800 
cm~^ in the northern sector of PA=169° (see the [N II] 
panel in Fig. 1). The match of A^e (forbidden lines) and 
A'^e(Ha) provides (rcspixDxei)outer shell ^^^10-12 arcsec kpc 
(rcspi= size of the outer shell in the radial direction), im- 
plying that (Dxei)outer shell —1.0 kpc, i.e. the same value 
derived for the main shell (Sect. 4.2). 

4.6. The caps 

In these medium-latitude condensations of the outer shell 
Te[0 III]~ Te[N II]~ 10 000 K, and 7Ve[S II] =^3000-3500 
cm~^; similar results are obtained by Balick et al. (1994), 
Kwitter & Henry (1998) and Rubin et al. (2002), whereas, 
according to Gongalves et al. (2003), Tc[0 III]~ rc[N 
II]~9300 10400 K, and iVe[S II]=4500-5000 cm^^. The 
large density reported by Gongalves et al. appears incom- 
patible with the weakness of the caps in the Ha panel 
of Fig. 1 (since, in first approximation, l{}la)(xNe ^). The 
same occurs when considering both the Ha and [O III] 
spectral images at PA=79° (Fig. 2, bottom panel). Our 
conclusion is that, due to projection effects, the line inten- 
sities provided by the low-resolution spectra of Gongalves 
et al. (2003; positions K2 and K3 in their Table 1) repre- 
sent a mix of both the caps and the underlying main shell, 
which is denser and at higher excitation (as confirmed by 
the analysis of our echellograms at PA=79° in "low spec- 
tral resolution" mode, i.e. integrated along the expansion 
velocity field). This is a clear demonstration of the severe 



limits implicit in the study of objects as complex as NGC 
7009 at unsuitable spectral resolution. On the contrary, 
the tomographic and 3-D methodologies, based on high- 
resolution spectra, reconstruct the true spatial structure 
of the gas and overcome any misleading camouflage due 
to projection. We stress that the matter distribution in 
Fig. 8 (bottom panel) enforces the indication outlined by 
the kinematics (Sect. 3.4): the low-excitation degree of the 
caps is likely due to the shadowing by the main shell, a 
triaxial ellipsoid of large ellipticity, quite dense along the 
radial directions containing the caps. 

4.7. The halo 

No To determination is possible for this weak sub- 
system, whose electron density peaks at the inner edge 
(~250-300 cm~^), and decreases outward. We assume 
Te(halo)~10000 K. 

5. The nebular distance, size, mass and age 

The determination of the distance, D(NGC 7009), through 
the interstellar absorption-distance relation for the field 
stars (with accurate photometry, spectral type and lu- 
minosity class) furnishes poor results, due to the spread 
of the data and the modest nebular extinction (Lutz 
1973, Gathier et al. 1986, Saurer 1995). A lower limit 
to the distance, D(NGC 7009)>1.0 kpc, is inferred from 

(Dxei) main shell - (D X ei)outer shell -1-0 kpC (SectS. 4.2 

and 4.5) and the assumption ei=l. 

Our nebula is listed in more than two dozen catalogues 
of "statistical" distances (since O'Dell 1962), providing 
the following mean values: 

<D>(Shklovsky)~1.3(±0.4) kpc 
<D>(ionized mass-radius relation) ~ 1.0 (±0.5) kpc 
<D>(surface bright .-radius relation)~1.0(±0.5) kpc 
<D>(other methods)~1.5(±0.6) kpc. 

NGC 7009 being optically thin in most directions, 
<D>(Shklovsky) is expected to be more or less reliable; 
the same is valid for <D> (ionized mass-radius relation) 
and <D> (surface brightness-radius relation), but at a 
lower degree of confidence, due to the difficulty in defining 
the radius of such an elongated object. 

Note that <D>(Shklovsky) disagrees with both "indi- 
vidual" distances reported in the literature: 

- Filler (1965) and Filler et al. (1966) measured the an- 

gular expansion of the inner ring, the outer ring and 
the ansae on different epoch (1899, 1908 and 1961) 
plates taken with the 36-inch Crossley reflector, ob- 
taining D(NGC 7009)~0.6 kpc; 

- Mendez et al. (1988) derived the mass, distance and lu- 

minosity of the central star by comparing the observed 
photospheric absorption-line profiles with the detailed 
profiles for a non-LTE model atmosphere. They give 
D(NGC 7009)~2.5 kpc (later revised to 2.1 kpc by 
Mendez et al. 1992). 
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Fig. 8. Radial distribution of the physical conditions (To and N^) in selected directions of NGC 7009: top panel= the 
zvpc at PA=169° (apparent minor axis), middle panel= the zvpc at PA=79° (apparent major axis), bottom panel= 
the de-projected direction passing through the central star and the caps in PA=94°. Left ordinate scale: Te[0 III] 
(continuous line) and To[N II] (squares). Right ordinate scale= Nc[S II] (circles), 7Vo[Ar IV] (dots) and Nc(Ra) for 
three representative values of eixrcspixD (dotted line= 5 arcsec kpc, short-dashed line= 10 arcsec kpc, and long-dashed 
line= 20 arcsec kpc). 



Let us enter in some detail. According to Liller (1965) 
and Liller et al. (1966), the nebular material moves out- 
ward at a velocity proportional to the distance from the 
central star. They adopt the value ^=0.7(±0.3)xl0-2 
arcsec yr^^ for the gas of the outer ring at 13 arcsec from 
the star; when combined with T4xp=20.5 km s~^ (Wilson 
1950), this yields D(NGC 7009)=0.6(lo:2) kpc, where: 



D(pc) 



0.21iyexp(km s-^) 
^(arcsec yr~^) 



(8) 



After 1966 this distance has been widely adopted, 
and NGC 7009 is listed among the calibrators in a 
number of statistical studies on PNe. Now the kine- 
matical results of Sect. 3 indicate that this must be 
revised. In fact, the expansion velocity of the gas in 
the outer shell at 13 arcsec from the star becomes 
Vexp -41 km s^i (from Eq. (4)), and the Liller (1965) 
and Liller et al. (1966) angular expansion, combined with 
Eq. (8), gives D(NGC 7009)=1.25(t°:!^) kpc. These au- 
thors measured also the angular expansion in the ansae, 
^=1.6(±0.3)xl0-2 arcsec yr"!; from Eqs. (3) and (8) 
and r(ansae)~25" we obtain Kxp (ansae) ~ 100 km s~^, 
and D(NGC 7009) = 1.30(t!^;^) kpc (quite uncertain, since 
in this case we combine the poorly known expansion veloc- 
ity of the gas with the angular expansion of the ionization 
front). 

The new values agree with <D>(Shklovsky), but are 
considerably lower than the "gravity distance" based on 
high-resolution spectra of the central star (Mendez et al. 
1988, 1992). The literature contains contrasting reports on 
the general validity of the "gravity distance" for the cen- 
tral stars of PNe (sec Pottasch 1996, Pottasch & Acker 
1998, Napiwotzki 1999, 2001). In the specific case of NGC 
7009, Mendez et al. (1988, 1992) obtained a quite large 
stellar mass (M*=0.70 and 0.66 Mq, respectively; well 
above the canonical mean value for PNe nuclei, 0.60 Mq; 
Blocker 1995 and references therein), whose fast evolution 
appears incompatible with the observations (see below). 
Note that the same contrasting result occurs (in the op- 
posite direction) when considering the distance suggested 
by Liller (1965) and Liller et al. (1966), i.e. 0.6 kpc: in this 
case the central star parameters correspond to an under- 
massive (M* <0.55 Mq), too slowly evolving post-AGB 
star. 

The definitive answer to the problem of the NGC 
7009 distance will come from the angular expansion mea- 
sured in first- and second-epoch HST imaging, combined 
with the detailed spatio-kinematical results of Sect. 3. 
At present forty-four broad-band and interference filter 



WFPC2 images (often multiple exposures) of the nebula 
are contained in the HST archives, secured at five epochs 
covering the interval 1996.32 to 2001.36. Two disturbing 
factors immediately appear: (i) the nebular size exceeds 
the field of the central (PC) camera in the WFPC2 mo- 
saic; (ii) the telescope pointing changes from run to run, 
i.e. the nebula is located in different parts of the detector. 
For example, in the 1996.32 frames only the western re- 
gion of NGC 7009 lies in the PC camera (pixel scale 0.046 
arcsec pix~^), whereas in the 2000.27 frames the whole 
nebula is in a lateral (WF) camera (pixel scale 0.100 arc- 
sec pix~^). 

All things considered (i.e. (1) the exciting star is the 
centring point of the first- and second-epoch observations, 
(2) the longer the temporal interval, the higher the angular 
expansion accuracy, and, mainly, (3) the frames "belong" 
to someone else), we limit the analysis to the F555 im- 
ages of the 1997.40 (program GO 6119; P.I. Howard Bond) 
and 2001.36 (program GO 8390; P.I. Arson Hajian) nms. 
Multiple exposures are co-added, corrected for optical dis- 
tortions, aligned and rotated using IRAF packages (Reed 
et al. 1999, Palcn et al. 2002). We measure the appar- 
ent shift of the main shell close to the major axis (in the 
W-SW sector, at a distance of 11.2 arcsec from the cen- 
tral star), obtaining ^=6.5(±2.0) x 10^'^ arcsec yr~^, and 
D(NGC 7009)=1450(;^^^) pc (from Eqs. (3) and (8)). 

This: 

(a) agrees with the results previously obtained from the 
Liller (1965) and Liller et al. (1966) ground-based an- 
gular expansion combined with the detailed spatio- 
kinematical model of Sect. 3; 

(b) overcomes the (legitimate) doubts and criticisms on 
the expansion parallaxes of PNe (Schonberner, private 
communication): this method compares matter veloc- 
ities with pattern velocities, and can seriously under- 
estimate the distance (up to a factor of two) when us- 
ing an inaccurate spatio-kinematical model, as clearly 
shown by the Liller (1965) and Liller et al. (1966) an- 
gular expansion of NGC 7009, simply combined with 
the expansion velocity given by Wilson (1950). 

In the following we will adopt D(NGC 7009)=1400 pc, cor- 
responding to a distance from the galactic plane l^j ~780 
pc. The resulting linear sizes of the different sub-systems 
are: r(main shell)~0.038x 0.073 pc, r(outer shell)~0.080 
pc, and r(ansae)~0.161 pc. The local filling factor in the 
nebula is ei ~0.7 (Sects. 4.2 and 4.5), and the ionized mass 
Mion— 0.18(±0.03) M0 (obtained in various ways: from the 
H/3 flux, the radio flux, and the observed A^e distribution; 
AUer 1984, Pottasch 1984, Osterbrock 1989, Turatto et al. 
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Table 4. NGC 7009: chronology desumed from the neb- 
ular kinematics 



Time 

(yr) 



Evolutionary phase 



logT, 
(K) 



to-7000(±1000) start of the superwind 3.7 

to-5000(±1000) end of the superwind 3.7 

to-2000(±500) start of the photo-ionization 4.4 

to-1000(±500) full ionization (main shell) 4.7 

to today 4.9 



2002). Mion is close to the total nebular mass, NGC 7009 
being optically thin in most directions. 

Concerning the kincmatical age, tkin oc R/V^xp, of the 
"large-scale" sub-systems, we have: 

- tkin(iiiain shell)~1650 yr, 

- tkin (outer shell) ~2200 yr, 

- tkin (halo) >8500 yr. 

We can link these observational times with the the- 
oretical ones coming from the most recent evolutionary 
models for PNe (Schonberncr ct al. 1997, Corradi et al. 
2000, Marigo et al. 2001, Villavcr ct al. 2002 and references 
therein); these predict the following chronological phases: 
AGB wind, superwind, transition time, photo-ionization, 
fast stellar wind, recombination and re-ionization (the last 
two phases mainly refer to a nebula powered by a massive, 
fast evolving star). Under the reasonable assumption that 
the observed halo of NGC 7009 corresponds to the AGB 
wind, we infer that tkin (halo) is an upper limit to the onset 
of the superwind, whereas tkin (outer shell) and tkin (main 
shell) are correlated to the ionization phase. 

When taking into account the probable dynamical his- 
tory of the gas, we obtain the rough chronological sequence 
of Table 4 (the last column contains an estimate of the stel- 
lar temperature, T*). These evolutionary times, derived 
from the nebular kinematics, must be checked with those 
coming from the observed characteristics of the central 
star, as performed in the next section. 

6. The central star 

The exciting star of NGC 7009 has mv= 12.80(±0.10) 
(Gathier & Pottasch 1988, Tylcnda ct al. 1991, CiarduUo 
et al. 1999, Rubin et al. 2002) and spectral type O (with 
H I, He I and He 11 absorption lines, Mendez et al. 1988). 

The stellar magnitude, combined with both the 
total H/3 nebular flux, log F(H/3)obs=-9.78(±0.03) 
mWxm-2 (Copetti 1990, Acker et al. 1992, Lame & 
Poggc 1996, this paper), and the flux ratio F(A4686 
A)/F(H/?)==0.18(±0.03) (Bohigas et al. 1994, Tylenda et 
al. 1994, Hyung & Aller 1995a, b, Lame & Pogge 1996, 
this paper), provides the following Zanstra temperatures: 
log(TzH I)= 4.83(±0.05) and log(TzHe II)=4.95(±0.05). 
The Zanstra discrepancy confirms that the nebula is 



optically thin to the H Lyman continuum, i.e. TzH 
I<T*~TzHe H, in agreement with the results by Gathier 

& Pottasch (1988), Patriarchi et al. (1989) and Mendez 
et al. (1992). Moreover, from D(NGC 7009)=1.4 kpc and 
c(H/3)=0.15 we obtain log L*/L0(TzH I)=3.50(±0.10) 
and log L*/L0(TzHe II)=3.70(±0.10) (using the bolomet- 
ric corrections by Schonberner 1981). 

In summary, the central star of NGC 7009 is hot 
(logT*~4.95) and luminous (log L*/Lq~3.70). Its po- 
sition in the log L*-log T« diagram, compared with 
the theoretical evolutionary tracks by Schonberner (1981, 
1983), Iben (1984), Wood & Faulkner (1986), Blocker 
& Schonberner (1990), Vassihadis & Wood (1994) and 
Blocker (1995), corresponds to a 0.60-0.61 M© post-AGB 
star in the hydrogen shell nuclear burning phase, whose 
evolutionary times fairly well agree with the kinematical 
results of Table 4. 

However, such a satisfactory picture, when analyzed in 
detail, presents a number of disturbing drifts: 

I) the uncertainties in the model evolution for post- 
AGB are quite large (e.g. semi-empirical AGB mass-loss 
law, treatment of the residual hydrogen-rich envelope, en- 
ergy release and neutrino losses in the fading to the white 
dwarf region); 

II) the estimated mass-loss rate for the central star 
of NGC 7009 spans the broad range 3.0x10"!° M© yr"! 
(Ccrruti-Sola & Perinotto 1985) to 1.0x10"*^ M© yr"! 
(Bombeck ct al. 1986); this is connected to the uncertain- 
ties in the stellar gravity, temperature and luminosity, and 
in the physical conditions, the ionic and chemical compo- 
sition of the wind; 

III) the diffuse X-ray emission of NGC 7009, arising in 
shocks at the interface between the fast stellar wind and 
the nebular shell, corresponds to a plasma at T~1.8xl0^' 
K (Guerrero et al. 2002), well below the post-shock tem- 
perature expected of the wind speed (Vcdgc— 2700 km s~^, 
Cerruti-Sola & Perinotto 1985). Different temperature- 
moderating mechanisms have been introduced to over- 
come this discrepancy, such as: coUimatcd oiitflows, adi- 
abatic cooling, direct mixing and/or heat conduction be- 
tween the shock-heated and the photo-ionized gas (for a 
review, see Soker & Kastncr 2003); 

IV) according to Guerrero et al. (2002), in NGC 7009 
the thermal pressure of the X-ray emitting gas is large 
enough to push outward the internal, low-density nebular 
regions, whereas no evidence of acceleration appears in 
our echellograms, for the high excitation "internal shell" 
either. As suggested in Sect. 3, a kinematical study of the 
(quite faint) [Ne V] emission at A3425 A is highly desir- 
able. The literature is desolately scanty in this field. We 
recall the classical paper of Wilson (1950), who measured 
V'exp[Ne V] in eight bright PNe, obtaining Vc^p[Ne Yjc^^O 
for five targets, and T4xp[Ne V]=(the lowest value among 
the observed ions) for the remaining ones. Even though 
Meaburn & Walsh (1980) reported the detection of a 800 
km s^! wide component to the [Ne V] A3425 A emission of 
NGC 6302 (the first direct evidence of radiatively ionized 
fast stellar wind), this result is questioned by Casassus 
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et al. (2000, 2002): echellograms in the [Ne V] A3425 A 
emission (R~80000) and the IR coronal Unes (R~20000) 
of very high-excitation ions (up to [Mg VIII], IP 225 cV) 
show that all lines of NGC 6302 are very narrow, with no 
evidence of: (a) broad wings from fast stellar winds, (b) a 
wind-blown central cavity filled with a hot rarefied plasma 
(also see Oliva et al. 1996); 

V) the radial density profile of the hydrodynamical 
PNe models (Icke et al. 1992, Mellema 1997, Schonberner 
et al. 1997, Dwarkadas & Balick 1998, Garcia-Segura et 
al. 1999, Prank 1999, Corradi et al. 2000, Blackman et 
al. 2001, Soker & Rappaport 2001, Marigo et al. 2001, 
Villaver et al. 2002) presents an internal, large, empty 
region (the hot bubble), whose bulldozer action on the 
nebular gas causes a very steep rise to the maximum den- 
sity. This disagrees with the quite smooth internal radial 
density distribution observed in NGC 7009 (this paper), 
NGC 1501 (Ragazzoni et al. 2001), NGC 6565 (Turatto 
et al. 2002), NGC 6818 (Benetti et al. 2003), and many 
other surveyed PNe. The odd thing is that, at present, 
the best match to the radial profiles of the observed PNe 
is provided by the hydrodynamical models neglecting (or 
limiting) the push effect by the fast wind on the innermost 
nebular regions; see, e. g., Bobrowsky & Zipoy (1989) and 
Frank et al. (1990), though: (a) Marten & Schonberner 
(1991) argue that Bobrowsky & Zipoy (1989) adopted un- 
realistic boundary conditions -no fossil AGB wind, no self- 
consistent central star evolution and inner rim pressure-, 
and Frank et al. (1990) did not consider the dynamical ef- 
fects of photo- ionization; (b) according to Mellema (1994), 
too low mass-loss rates (or velocities) in the fast post- AGB 
wind can cause the collapse of the nebula when ionization 
sets in (but see also Chevalier 1997). 

The matter gets further complicated when introducing 
the "small-scale" sub-systems of NGC 7009 (the caps, the 
streams and the ansae). We will search for new input in 
the observed radial ionization structure and the photo- 
ionization model. 

7. The mean chemical and radial ionic 

abundances, and the photo-ionization model 

The mean ionic abundances of NGC 7009 arc obtained 
from the line fluxes integrated over the spatial profile and 
the expansion velocity field, according to the critical anal- 
yses by Alexander & Balick (1997) and Pcrinotto et al. 
(1998). Once the ionization correcting factors for the un- 
observed ionic stages have been considered (derived both 
empirically. Barker 1983, 1986, and from interpolation of 
theoretical nebular models. Shields et al. 1981, AUer & 
Czyzak 1983, AUer 1984, Osterbrock 1989), we find the 
chemical abundances listed on the bottom line of Table 5 
(to be compared with the estimates reported in the liter- 
ature, also contained in the table). They are mean values 
representative of the main, bright nebula, and have been 
used (in Sect. 4) to derive F(Ha)zvpc and Nc{lla). 

The detailed radial ionization structure of NGC 7009 
is given by the zvpc of the different lines. Since the large 



Ha broadening limits the accuracy of the F (Ha) zvpc dis- 
tribution, we adopt A5007 A as reference emission, thus 
inferring the ionization structure relative to 0++ through: 



F{X{X+''))^ 



0++ i^(A5007A), 



(9) 



where a=b for the forbidden lines, and a=b-|-l for the 
recombination lines (for details, see Benetti et al. 2003). 

The -^^r^ profiles in the zvpc at PA=169° (apparent 
minor axis), the zvpc at PA=79° (apparent major axis), 
and the de-projected direction containing the central star 
and the caps in PA=94° (see Sect. 4) arc presented in Fig. 
9. They confirm and/or enforce many of the observational 
results outlined in the previous sections. In particular: 

- along the minor axis (top panel) the nebula is optically 

thin to the UV stellar flux. Note, however, the prc^sence 
of a dense, quasi-thick condensation in the northern 
sector of the outer shell (clearly visible in Fig. 1); 

- along the major axis (middle panel) the internal shell, 

the main shell and the streams are high-excitation re- 
gions. The ionization quickly decreases in the ansae, 
whose outer edges are characterized by strong [O I] 
emission, i. e. they are neutral (or almost neutral) 
zones dominated by the resonant charge-exchange re- 
action 0+ + H0t^O° + H+ (Williams 1973); 

- the same ionization drop occurs in the caps (bottom 

panel), which are shaded by the dense main shell. 



Note that in Fig. 9 both the and radial profiles 

tend to increase inward, suggesting an overabundance of 
He in the faint, innermost regions of NGC 7009. Although 
weakened by the low line fluxes and the He II A4686 A fine- 
structure broadening, this agrees (at least qualitatively) 
with the different nebular morphology at A4686 A and 
Ha (Fig. 1). Let us enter in more detail, profiting of the 
superior spatial resolution of HST. 

Fig. 1 presents a large, well-defined stratification of 
the radiation: He+"'" marks the highest-excitation layers, 
represented by the central hollow and the inner part of the 
main shell, whereas H"*" is emitted by the whole nebula, i.e. 
the central hollow, the entire main shell, the outer shell 

and the halo. In order to obtain the intrinsic ^^^fff^^^ 

I{Ha) 

ratio in the internal, high-excitation layers of NGC 7009, 
the Ha contamination by the halo, the outer shell and the 
external part of the main shell must be removed in Fig. 1. 
This can be done in two ways: 

(a) from simple geometrical considerations, assuming 
spherical symmetry for the different emitting strata; 

(b) using the F502 image as a marker of the He"'""'"- 
poor layers, since IP(0++)=54.9 eV~IP(Hc+)=54.4 eV; 
i.e. the [O III] A5007 A line identifies the regions where 
He is single-ionized. 

Both these methods yield ^^^^^^^ >1.5 in the inner- 
most part of the He"*""*" shell and the central hollow of Fig. 
1. Now, the Te dependence of these two recombination 

lines is weak, and comparable. For a given ^^ j^^^^^ value, 
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Reference 


Position 


He 


C 


N 


O 


Ne 


S 


01 


Ar 


Liu et al. 1995 


A 


11.04^*) 


8.96^*^ 


8.71^*' 


9.28^*' 










Hyung & AUer 1995a,b 


B 


10.99 


8.43 


8.26 


8.73 


8.13 


7.14 


5.26 


6.40 




A 


11.07 


8.12 


8.06 


8.58 


8.12 


6.88 


5.22 


6.42 




G 


11.04 


8.66 


8.46 


8.76 


7.91 


- 


- 


- 


Kwitter & Henry 1998 


A 


11.08 


8.70 


8.55 


8.73 


7.90 


- 


- 


- 




E 


11.08 


8.61 


8.48 


8.76 


8.09 










H 


11.04 


8.62 


8.52 


8.74 


7.97 








Luo et al. 2001 


A 


- 


- 


- 


- 


8.24 


- 


- 


- 




A 




- 


- 


- 


8.84(*) 


- 


- 


- 




B 


11.03-11.07 




7.84-8.26 


8.65-8.68 


8.05 


6.69-6.78 






Gongalves et al. 2003 


C+H 


11.02-11.04 




8.26-8.38 


8.67-8.79 


8.05-8.08 


6.91-7.20 








S 


11.03-11.09 




7.79-7.93 


8.48-8.81 


7.98-8.05 


6.64-6.80 








D+E 


10.98-11.01 




8.40-8.58 


8.65-8.76 


8.05-8.11 


6.97-7.14 






This paper 


F 


11.05 


8.2: 


8.21 


8.70 


8.15 


7.15 


5.3: 


6.40 



(*) from optical recombination lines A=main shell (minor axis) D=West ansa G=outer shell 

B=main shell (major axis) E=East ansa H=East cap 

C=West cap F=overall nebula S=streams 



Fig. 9. The radial ionization structure of NGC 7009 in the zvpc at PA=169° (apparent minor axis; top panel), 
the zvpc at PA=79° (apparent major axis; middle panel), and the dc-projected direction passing through the central 
star and the caps in PA=94° (bottom panel). Same orientation as Fig. 8. 



an increase of Te implies a lower increase of ^4. (much 
lower for rc>15 000 K), and an almost parallel cmissivity 
decline (for example, when doubles, going from 10 000 
to 20000 K, increases by 12%, and the emissivities 

decrease by 53% for A4686 A, and 47% for Ha; AUer 1984, 
Pottasch 1984, Ostcrbrock 1989). Thus, in first approxi- 
mation we can put ^^X^ — 0.23 x ^^jjjf^y^- Moreover, in 
the internal, high-excitation regions — ^f- 

We obtain '^+^ >0.30 in the internal part of the 
He++ shell and the central hollow of NGC 7009 (well 
above the mean helium abundance listed in Table 5). 
Similar (but less precise) results come from our cchel- 
lograms, when taking into account the A4686 A and Ha 
intensity profiles integrated over the whole expansion ve- 
locity field. All this, on the one hand supports the suspi- 
cion raised by the radial ionization structure, but on the 
other hand is far from conclusive, due to the large Ha 
contamination by the halo, the outer shell and the exter- 
nal part of the main shell. The literature is scarce and 
controversial: concerning NGC 7009, echellograms of the 



exciting star by Mendez et al. (1988) provide ^=0.05 at 
the stellar surface, and Guerrero et al. (2002) fit the X- 
ray spectrum of the central cavity using normal nebular 
abundances, whereas the extended X-ray spectrum of: 

- NGC 6543 needs an overabundance of He by a factor 
of 60 with respect to the solar value (Chu et al. 2001; this 
result is questioned by Maness & Vrtilek 2003); 

- BD-|-30°3639 presents a significant continuum emis- 
sion, suggestive of greatly enhanced He abundances 
(Kastner et al. 2000); 

- NGC 7027 is overabundant in He, C, N, Mg and Si 
(Kastner et al. 2001). 

Moreover: 

- HST NICMOS imaging of NGC 5315 in different 
emission lines allowed Pottasch et al. (2002) to determine 
the helium abundance across the nebula: it is higher on the 
inside than on the outside, but not uniformly distributed 
(i.e. the He-rich gas appears in two opposite directions); 

- both the kinematics and the radial ionization profile 
of NGC 40 suggest that the fast, H-depleted photospheric 
material ejected by the WC8 central star is gradually mod- 
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ifying the chemical composition of the innermost nebular 
regions (Sabbadin et al. 2000a). 

7.1. Photo-ionization model 

In this section the photo-ionization code CLOUDY 
(Ferland et al. 1998) is applied to a PN characterized by 
the same distance, gas distribution, mean chemical com- 
position and exciting star parameters of NGC 7009. As 
usual (Sects. 4 and 7), we take into account different ra- 
dial directions: 

- the zvpc in the S-SE sector at PA=169° for the main 

shell, the outer shell and the halo along the apparent 
minor axis, 

- the zvpc in the E-NE sector at PA=79° for the internal 

shell, the main shell, the streams and the ansae along 
the (almost true) major axis, 

- the de-projected direction containing the central star 

and the bright cap in the West sector at PA=94°. 

The input parameters of the model nebula are given in 
Table 6, whereas Fig. 10 presents the absolute radial flux 
distribution of the main emissions in the model nebula 
and the real nebula. We infer the following crucial points: 

(I) all the spectral features observed in the sub-systems 
of NGC 7009 are fairly well explained by photo- 
ionization; an even better fit is obtained simply lower- 
ing ei(model) (0.68 instead of 0.70), and/or T*(model) 
(85 000 K instead of 90 000 K). 

This: (a) agrees with the results by Balick et al. (1994, 
1998), Hyung & Aller (1995a, b) and Gongalves et al. 
(2003), indicating that the nebula is mainly radiativcly 
excited; (b) doesn't exclude that some shocked regions 
can be present in the ansae (as expected from the su- 
personic motion): following Dopita (1997), in the case 
of a shock illuminated by an external photo-ionizing 
source, the energetic flux of the ionizing photons over- 
whelms the mechanical energy flux through the shock; 

(II) no N-enrichment need be postulated in the caps and 
the ansae to account for the local enhancement of the 
low excitation emissions. This feature is a mere con- 
sequence of the ionization drop: when increasing the 
gas distance from the central star, the ionization equi- 
librium moves towards lower states since the recombi- 
nation rate is almost unchanged, whereas ionization is 
less frequent (more details are in Alexander & Balick 
1997 and Gruenwald & Viegas 1998). At aU radial di- 
rections of NGC 7009 (excluding the major axis) the 
main shell is a dense, optically thin layer degrading 
most of the UV stellar radiation, so that the mean 
latitude condensations in the outer shell (i.e. the caps) 
have to be optically thick. The main shell is weak along 
the major axis, and the UV stellar flux is here degraded 
in extended regions of moderate density (the streams), 
whose external, slightly denser edges (the ansae) yield 
the rapid plasma recombination. 



Table 6. Input parameters for the model nebula 
(CLOUDY) 



Radial density profile 


Fig. 8 (cfr. Sect. 4) 


Chemical abundances: 




C, F, CI, K, Ca 


Hyung & Aller (1995a,b) 


He N O Ne S Ar 




irf" n PT* plpTnPTi't'Q 


PN (^CT.OTinV default "1 


Dust 


PN (CLOUDY default) 


Local filling factor 


0.7 




blackbody distribution 


Exciting star 


T»=90000 K 




log L,/L0= 3.70 


Distance 


1.4 kpc 



Fig. 11. Stereoscopic structure of NGC 7009 for a rota- 
tion around the North-South axis centered on the exciting 
star. Opaque reconstruction (at high and low flux cuts) for 
A4686 A of He II, as seen from seven directions separated 
by 15°. The line of view is given by {0,tp), where 6 is the 
zenith angle and ijj the azimuthal angle. Each horizontal 
pair represents a "direct" stereoscopic pair (i. e. parallel 
eyes), and the whole figure provides six 3-D views of the 
nebula (at two flux cuts) in as many directions, covering 
a right angle (for details, see Ragazzoni et al. 2001). The 
(0,0) images represent the rebuilt-nebula seen from the 
Earth (North is up and East to the left). 



8. The 3-D morpho-kinematical structure 

The 3-D reconstruction method for expanding nebulae 
was introduced by Sabbadin et al. (1985, 1987), and re- 
fined by Sabbadin et al. (2000a, b). In the case of the 
Saturn Nebula we select A4686 A of He II, A5007 A of [O 
HI] and A6584 A of [N II] as markers of the high, mean 
and low- ionization regions, respectively. 

The spectral images of the forbidden lines are de- 
convolved for seeing, spectral resolution and thermal mo- 
tions, whereas also fine-structure is taken into account for 
A4686 A of He II. The emissions of the main shell (+ the 
streams -|- the ansae) are de-projected through Eq. (3), 
and those of the outer shell (-1- the caps) through Eq. (4), 
and later assembled by means of the 3-D rendering proce- 
dure described by Ragazzoni et al. (2001). 

For reasons of space, in this paper we extract, present 
and discuss a limited number of frames, corresponding 
to a partial rotation around the North-South axis cen- 
tered on the exciting star (i.e. almost perpendicular to 
the major axis). They provide a good picture of all as- 
pects of the nebular complexity outlined in the previous 
sections. The complete series of nebular movies is avail- 
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Fig. 10. Physical conditions and radial ionization structure in selected directions of the model nebula (CLOUDY), 
and the true nebula. Left column: zvpc in the S-SE sector at PA=169° (apparent minor axis); central column: zvpc 
in the E-NE sector at PA=79° (apparent major axis); right column: along the de-projected direction containing the 
central star and the West cap in PA=94°. The top and the middle panels concern the model nebula, and provide the 
physical conditions and the absolute fluxes in the main emissions, respectively; the bottom panels present the absolute 
radial flux distribution observed in the true nebula. 



able at http: / /web. pd.astro.it / sabbadin' this WEB 
site dedicated to the spatial structure of the PNc must 
be considered as an integral (and important) part of the 
paper. 

Figs. If to f3 contain the opaque reconstruction of 
NGC 7009 in He II, [O III] and [N II] (at different flux 
cuts) for a rotation of 90° through the first Euler angle. 
The (0,0) images correspond to the Earth-nebula direc- 
tion (North is up and East to the left). 

The highest-excitation regions (Fig. If) form an in- 
homogeneous ellipsoid, bright in two equatorial caps and 
open-ended along the major axis; at lower A4686 A fluxes 
they merge into a closed structure with polar, diffuse lay- 
ers. 

The opaque reconstruction in [O III] (Fig. 12) is given 
for three absolute flux cuts: log E(A5007 A)= -16.82, - 
17.54 and -18.46 erg s~^ cm~'^ (high, mean and low cuts, 
respectively). Since O/H=5.0xl0-''~O++/H+ (Table 5 
in Sect. 7) and Nc ~1.15xN(H+), they correspond to 
A'^e(higli cut)~5300 cm"^, A^c(niean cut)~2200 cm~^, and 
A^e(low cut)~700 cm-3 (assuming Tc^lO 000 K and ei=l). 

The different sub-systems are clearly identifiable in 
these frames: the brightest regions of the main shell appear 
at high cut, the entire main shell, the caps and the densest 
parts of the outer shell at mean cut, and the overall nebu- 
lar structure, including the streams and the ansae, at low 
cut. Moreover, the mean cut frames in Fig. 12 highlight 
the presence of a further (the eighth!) sub-system: an ir- 
regular, equatorial pseudo-ring within the outer shell. The 
newly discovered ring is barely visible (almost edge-on) in 
the Ha and [N II] panels of Fig. 1 ; its most conspicuous 
feature is the quasi-thick knot in the northern sector at 
PA=169° (Sects. 4.5 and 7). 

In Fig. 14 the axes of the bright low-excitation sub- 
systems of NGC 7009, i.e. the ansae and the caps, are mis- 
aligned by 30°. At lower [N II] fluxes, the external parts of 
the main shell and the (just identified) equatorial pseudo- 
ring within the outer shell show up. Please note the spatial 
correlation between the caps (Fig. 13, high cut) and the 
densest regions of the main shell (Fig. 12, high cut), con- 
firming that the former do suffer strong shadowing by the 
latter (Sects. 3.4, 4, 7 and 7.1). 

The assembled, multi-color projection of NGC 7009 for 
a rotation around the N-S direction (almost perpendicular 
to the major axis) is presented in Fig. 14, providing a rep- 
resentative sample of the nebular appearance when chang- 
ing the line of view. The upper-right panel, (0,0), corre- 
sponds to NGC 7009 as seen from the Earth (North is up 
and East to the left), to be compared with Fig. 1 and the 



HST multi-color images by Balick ct al. (1998) and Hajian 
& Tcrzian at http://ad.usno.navy.mil/pne/gallcry.html 

Note the florilegium of silhouettes present in Fig. 14: 
actually, the Saturn Nebula covers most of the morpholog- 
ical classes, i.e. roundish to elliptical, to bipolar (according 
to Greig 1972, Stanghellini et al. 1993, Corradi & Schwarz 
1995, and Corny et al. 1997). A similar variety of forms 
occurs for the other PNe so far analyzed with the 3-D pro- 
cedure (i.e. NGC 6565 and NGC 6818; Turatto et al. 2002 
and Benetti et al. 2003, respectively; see the movies at 
http://web.pd.astro.it/sabbadin). Such a dominant role 
of the perspective questions the validity of the statistical 
analyses based on a mere morphological classification of 
PNe, and stresses the need for a new code founded on 
"intrinsic, physical" properties of the star+nebula system 
(instead of "rough, apparent" characteristics). 



9. Discussion 

We have successfully applied the 3-D method (valid for 
all types of extended, expanding nebulae; Benetti et al. 
2003 and references therein) to the complex PN NGC 
7009, covered at high spatial and spectral resolutions with 
ESQ NTT-hEMMI. The Saturn Nebula hes at a distance 
of 1.4 kpc (age ~6000 yr, ionized mass ~0.18 M©), and 
consists of four "large-scale" sub-systems (the internal 
shell, the main shell, the outer shell and the halo) and 
as many "small-scale" ones (the streams, the caps, the 
ansae and the equatorial pseudo-ring), all interconnected 
and characterized by different morphology, excitation de- 
gree, physical conditions and kinematics. The four "large- 
scale" sub-systems and the streams are medium-to-high- 
excitation, optically thin regions, the caps and the ansae 
are thick, low-to-very low-excitation layers, whereas the 
equatorial pseudo-ring within the outer shell is quasi-thick 
(medium-to- low-excitation) . The internal shell, the main 
shell, the streams and the ansae follow the expansion law 
V^xp=4.0xR" km s^^, the outer shell, the caps and the 
pseudo-ring move at Vcxp— 3.15xR"km s^^, the halo ex- 
pands at Vcxp — fO km s^^. 

The accurate definition of the 3-D ionization structure 
of NGC 7009 allows us to combine the radial distribution 
of the physical conditions and the line fluxes in the dif- 
ferent sub-systems with the model proflles given by the 
photo-ionization code CLOUDY. We infer that all the ob- 
served features are explainable in terms of ionization by 
the central star, which is a hot (logT,~4.95) and lumi- 
nous (log L,/L0~3.7O) 0.60-0.61 M© post-AGB star in 
the hydrogen-shell nuclear burning phase. 
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Fig. 12. Opaque reconstruction of NGC 7009 at A5007 A of [OIII] (higli, mean and low flux cuts). Same scale as Fig. 
11. 

Fig. 13. Opaque reconstruction of NGC 7009 at A6584 A of [Nil] (high and low flux cuts). Same scale as Figs. 11 and 
12. 



At this point we can compare the overall characteris- 
tics of NGC 7009 (star + nebula) with the two other PNe 
so far analyzed in detail with the 3-D methodology: 

NGC 6565 (Turatto et al. 2002) is a young (2300 yr), op- 
tically thick ellipsoid embedded in a large cocoon of 
neutral, diisty gas. It is in a deep recombination phase 
(started about 400 yr ago), caused by the luminos- 
ity drop of the massive powering star (M*>O.65M0), 
which has almost reached the white dwarf region (log 
L,/Lq~2.0, log T«~5.08); 

NGC 6818 (Benetti et al. 2003) is a young (3500 years), 
optically thin (quasi-thin in some directions) double 
shell seen almost equator-on. It is at the start of the 
recombination phase following the luminosity decline 
of the 0.625 Mq central star, which has recently ex- 
hausted the hydrogen shell nuclear burning and is 
rapidly moving toward the white dwarf domain (log 
-5.22 K; log U/Lq -3.1). 

NGC 7009 is found to be the oldest, but also the least 
evolved, member of the trio, confirming that the stellar 
mass is the main driver of PNe evolution (Blocker 1995 
and references therein). 

When searching for more examples of NGC 7009-like 
PNe in the imagery catalogues (Acker et al. 1992, Schwarz 
et al. 1992, Manchado ct al. 1996, Gorny et al. 1999, and 
the HST archives), perspective effects play a leading role 
(as indicated by Fig. 14). We suspect this is also the case 
for NGC 2392, the famous "Eskimo Nebula". According 
to Reay et al. (1983), O'DeU & Ball (1985), O'DcU ct al. 
(1990) and Phillips & Cuesta (1999), it consists of an in- 
ner, bright, fast expanding (— 93 km s^^) prolate spheroid 
(axial ratio 2:1) much denser in the equatorial regions, 
seen almost pole on. A high-velocity stream (— 190 km 
s""'^) is formed at the open ends of the inner shell, proba- 
bly connected to the accelerating effects by the fast stellar 
wind. NGC 2392 is embedded in a slow expanding (~ 16 
km s~^) spheroid containing a fringe of fast (— 75 km 
s~^), low-excitation knots and filaments. Moreover, a col- 
limated, jet-like bipolar flow with T^xp —200 km s^^ has 
been detected by Gieseking et al. (1985) at PA=70°. The 
exciting star of the Eskimo Nebula has TzH 1=32 700 K, 
TzHe 11=73 700 K, log L*/Lq=3.99, and a mass-loss rate 
below 8.0x10-1° M© yr"!, with Vedge-400-600 km s'^ 
(assuming a distance of 1.5 kpc; Tinkler & Lamers 2002). 

All this suggests the structural similarity of NGC 2392 
and NGC 7009: the former is projected almost pole-on 
(see Fig. 9 in O'Defl & Ball 1985), the latter almost 
equator-on. Moreover, NGC 2392 is in a slightly earlier 
evolutionary phase, i.e. the main shell is fully ionized, but 
the outer shell is still thick (or quasi-thick) to the UV 



stellar flux, as supported by a preliminary analysis of the 
TNG-hSARG (eight PA, R-115 000) and Asiago+Echelle 
(twelve PA, R-25 000) spectra of the Eskimo Nebula. 

Further examples of NGC 7009-like PNe are IC 4634, 
J 320, NGC 2371-2 and HB 4. 

Now we are ready to combine the large amount of kine- 
matical and physical information stored in the previous 
sections with the theoretical predictions (Blocker 1995, 
Mellema 1997, Capriotti 1998, Garcia-Segura et al. 1999, 
Marigo et al. 2001, Villaver et al. 2002 and references 
therein), and sketch the evolution of an NGC 7009-like 
PN in a series of snapshots. 

The flrst one, at time to-5000 yr, corresponds to the 
end of the superwind (this enhanced mass-loss phase 
lasts 1000-2000 yr). The central star parameters are: 
logT*-3.70 and logL^/L® -3.75. The newly born PN con- 
sists of about 0.15 Mq of slowly expanding (Voxp=10-20 
km s^^) neutral gas (-1- molecules -|- dust), whose struc- 
ture is the result of two different superwind regimes: the 
older, almost spherical and isotropic, has been followed 
by a stronger, polar-depleted mass-loss. Moreover, a large 
cocoon of AGB wind embeds the nebula. 

In the period to-5000 yr (logT*~3.70 and logL^/Lg 
-3.75) to to-2000 yr (logT*-4.40 and logL*/L0~3.75), 
corresponding to the transition phase, the UV stellar flux 
is still negligible, and the neutral proto PN expands in a 
homologous fashion (Voxp=10 20 km s^^) maintaining the 
large pole-on to pole-off density contrast, also sustained by 
the moderate- velocity stellar wind, which prevents the gas 
backflow. According to Canto et al. (1988) and Frank et 
al. (1996), the stellar wind can, eventually, be deflected 
poleward by interaction with the ellipsoidal nebula. 

Around to-2000 yr (logT*-4.40 and logL,/L0-3.75) 
photo-ionization starts, increasing the temperature (i.e. 
the kinetic energy) of the innermost nebular layers, which 
expand untill they reach pressure balance with the sur- 
roundings: inward expansion is prevented by the thermal 
pressure of the faster and faster stellar wind, whereas the 
external, adjacent, neutral regions are pushed outward, 
inducing fragmentation (i.e. Rayleigh-Taylor instabilities; 
Kahn & Breitschwerdt 1990, Breitschwerdt & Kahn 1990). 
Since the propagation velocity of the Stromgren radius de- 
pends on both the UV stellar flux and the density of the 
neutral gas, in all radial directions of NGC 7009 (except 
the major axis) the ionization front is trapped inside the 
dense main shell, whereas along the major axis it soon 
penetrates the outer shell (Mellema 1995). 

At to-1000 yr (logT*-4.70 and logL*/LQ-3.75) the 
ionization front includes the whole main shell, and in all 
directions (excluding the major axis) the gas pressure acts 
on the outer shell, which is still neutral. Along the major 
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Fig. 14. Multi-color appearance of NGC 7009 (blue=He II, green=[0 III], red=[N II]) for a rotation through the N-S 
axis centered on the exciting star. The upper-right panel, (0,0), corresponds to the re-built nebula as seen from the 
Earth (North is up and East to the left). Recall that projection(^, V')=projection(^ ± 180°, ± 180°). 



axis the ionization front has reached the external edge of 
the outer shell by now, further accelerating into the halo 
(shocks?). Note that in this phase the radiative cooling 
is inadequate, due to the increased velocity of the stel- 
lar wind, and an adiabatic shock forms at the interface 
with the nebular gas, squeezing the innermost layers of 
the main shell. Due to the compression of the main shell, 
the recombination rate increases, thus retarding the time 
at which the main shell is fully ionized (Breitschwerdt & 
Kahn 1990). 

In the last evolutionary snapshot, referring to to 
(present time; logT«~4.95 and logL*/L0~3.7O), the whole 
nebula is optically thin to the UV stellar radiation, with 

the exception of: 

- the supersonic ionization fronts along the major axis, 

located well within the halo (i.e. the ansae), 

- two strings of symmetric, mean latitude, thick conden- 

sations still present in the outer shell (the caps). 

What next? At present the post-AGB star has reached 
the flat top of the UV flux (Blocker 1995 and references 
therein), and in the next millennium we do not expect 
great changes in the nebular structure, since the gas di- 
lution due to expansion will balance the smooth decline 
of the ionizing photons. Later on, the stellar evolution 
presents a quick drop in luminosity (at the end of the 
hydrogen-shell nuclear burning), and the recombination 
processes will prevail on photo-ionization within the neb- 
ula. The recombination time scale being proportional to 
^ (Aller 1984, Pottasch 1984, Osterbrock 1989), the caps 
and the ansae will be the first to rccombine (roughly at 
to+1100 yr; the ansae will showily contract), followed by 
the outer shell and the external layers of the main shell (at 
to+1200 yr), whereas the low-density halo will maintain 
a high-excitation degree for a long time. Beginning from 
to+2000 yr (logT^^S.l and logL,/Lo~2.5), the slower 
and slower stellar decline, combined with the gas dilution 
due to expansion, will bring on the gradual re-ionization 
of the nebula (of course, starting with the outer layers of 
the main shell). 

Summing up: 

- photo-ionization is the dominant process of the NGC 

7009 shaping (but not the exclusive process, since the 
fast stellar wind models the innermost nebular regions; 
see Capriotti 1973, 1998). In particular, the peculiar 
excitation, structure and kinematics of the streams 
and the ansae recall the "champagne-phase" model 
(Tenorio-Tagle 1979) for H II regions with an outward 
decreasing density distribution: when the density gra- 
dient exceeds some critical value, the pressure gradi- 
ents produce a weak R-type shock front that accel- 
erates the gas to supersonic velocities (Franco et al. 



1990, Garcia-Segura & Franco 1996, Shu et al. 2002); 
according to Spitzer (1968), the speed of a weak R- 
type front is nearly an order of magnitude larger than 
the velocity of a Ttypc one (which is about the sound 
speed in the ionized region, i.e. ~10 km s~^); 

- the outer shell and the caps of NGC 7009 are kinemat- 

ically coeval (the latter simply representing local con- 
densations within the former; Sects. 3.4 and 8); this is 
correlated to the last ionization phase of the nebula; 

- the three further co-moving sub-systems (i.e. the main 

shell, the streams and the ansae) can be considered 
kinematically coeval too, being associated to the start 
of the ionization process; 

- the agent(s) and mechanism(s) producing the two-phase 

superwind ejection (first isotropic, then deficient along 
the major axis) are still nebulous: jets in the pre- and 
proto-PN phases, magnetic fields, binary central star, 
accretion disc around a close companion (for details, 
see Sahai & Traugcr 1998, Sokcr & Rappaport 2000, 
Gardiner & Frank 2001 and references therein); 

- detailed 1-D and 2-D hydrodynamical simulation are 

needed (and highly desired) to improve/reject the 
rough picture here outlined. 

Our simple evolutionary model (nearly a challenge), 
on the one hand interprets fairly well the observational 
characteristics of NGC 7009, on the other hand does not 
pretend to be the panacea for all microstructures of PNe. 
Gongalves et al. (2001) list 50 objects with LIS (low ion- 
ization structures, including jets, jet-like systems, sym- 
metrical and non-symmetrical knots), analyze the mor- 
phological and kinematical properties, and compare them 
with the current theoretical models. An excellent review 
on FLIERs, jets and BRETs (bipolar rotating episodic 
jets) is provided by Balick & Frank (2002, and references 
therein). Moreover, Corradi et al. (2000) and Bcnetti et al. 
(2003) demonstrate that the low ionization structures ex- 
hibited by NGC 2438 and NGC 6818 can be explained in 
terms of a nebula in the recombination phase (caused by 
the luminosity decline of the hot central star, which has ex- 
hausted the hydrogen shell nuclear burning and is moving 
toward the white dwarfs domain) . We believe this is a not 
uncommon situation among PNe with LIS (cf. NGC 2440, 
NGC 2452, NGC 7293 and NGC 7354; Tylenda 1986). 

Let us focus on the ansae and the caps of NGC 7009. 
Strictly speaking, neither of them can be considered as 
genuine FLIERs, i.e. localized regions, whose characteris- 
tics (low excitation, N overabundance and peculiar motion 
with respect to the expansion law of the nebula) are sug- 
gestive of an exotic origin and/or evolution. In fact: 

- the caps are co-moving with their surroundings (i.e. the 

outer shell), 
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- the ansae follow the same kinematics of the main shell 

and the streams, 

- the low excitation degree of both sub-systems can be ex- 

plained in terms of "normal" nebular gas photo-ionized 
by the UV flux of the central star. 

Only in a broad sense are the ansae of NGC 7009 FLIERs, 
because of the peculiar motion with respect to the halo. In 
this case the streams become the prototype of a brand-new 
class of small-scale structures: the (sic!) LIFERs (large 
ionization, fast emitting regions). 

Joking aside, the LIS in PNe remain a tangled skein. 
Given the variety of phenomenologies, and of physical and 
evolutional processes entering into the picture, we con- 
clude that only a detailed analysis of a representative sam- 
ple of targets can provide reliable and deep insights. To 
this end we have collected high-quality echellograms for a 
dozen objects with LIS. 

Let us close the long journey dedicated to the Saturn 
Nebula listing the many problems left unsolved. 

Very likely the ansae expand supersonically, implying 
the presence of shocks (Sects. 3.1 and 7.1), but: (a) the me- 
chanical energy flux through the shock is overwhelmed by 
the energetic flux of the photo-ionization (Dopita 1997); 
(b) the motion is almost tangential (i.e. in the plane of 
the sky) , and the spatio-kinematical analysis based on the 
(small) radial velocity component inaccurate (Sects. 3.1 
and 3.4). The same is valid for the streams. Deeper obser- 
vations at even higher spatial and spectral resolutions (and 
in a wider spectral range) will shed light on this and on 
three more fields (the first two only grazed in the paper, 
the third not even touched): (i) the traverse, turbulent 
motions (i.e. the FWHM) in the streams and the ansae 
(Sects. 2 and 3.1), (ii) the spatio-kinematics of both the 
innermost, very high excitation layers (Sects. 3.1 and 6), 
and the extended, faint emissions surrounding the streams 
and the ansae (Sect. 3.3), and (iii) the discrepancy in the 
chemical composition of the gas obtained from optical re- 
combination and coUisionally excited lines (Hyung & AUer 
1995b, Liu et al. 1995, Luo et al. 2001). 

We infer (Sect. 2) that: (a) the baricentric radial veloc- 
ity is the same in the outer shell, the caps, the ansae and 
the halo of NGC 7009, but these sub-systems are blue- 
shifted (by 3(±1) km s""'^) with respect to the main shell; 
(b) the motion of the streams is hybrid; (c) the tip of both 
ansae is blue-ward. Can, all this, be connected to a binary 
nature of the central star (as suggested by the infrared 
photometry of Bentley 1989)? Does the different spatial 
orientation of the ansae and the caps (Sect. 8) support 
the precession hypothesis for the central star, eventually 
driven by a close companion (Bohigas et at. 1994)? 

Moreover: how can we reconcile the observed and the 
theoretical radial density profiles in PNe (Sects. 4 and 6)? 
Does the He overabundance in the innermost regions of 
NGC 7009 (Sect. 7) extend to other elements (and to other 
nebulae), implying a revision of the stellar mass-loss rates 
(so far obtained assuming nebular or solar abundances in 
the fast wind)? Do the detailed hydrodynamical simula- 



tions confirm/discard the naive evolutionary model here 
proposed (Sect. 9) for the Saturn Nebula? What are the 
effects of a clumpy structure on the nebular evolution? 
And of magnetic fields? 

To answer such a crescendo of questions, thorough in- 
vestigations are required, both theoretical and observa- 
tional, in a wide series of specific fields. Thus, mindful 
of Apelles's warning: "Cobbler, stick to thy last!", we can 
only pass the baton to the specialists, and stop. 
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